o 

(N 



=3 



> 



O 



X 



a 



A PROBABILISTIC VIEW ON THE GENERAL RELATIVISTIC 
BOLTZMANN EQUATION 

I. BAILLEUL 1 AND F. DEBBASCH 2 



Abstract. A new probalistic approach to general relativistic kinetic theory 
is proposed. The general relativistic Boltzmann equation is linked to a new 
Markov process in a completely intrinsic way. This treatment is then used to 
prove the causal character of the relativistic Boltzmann model. 



r*j ' 1. Introduction 

di- 

I Relativistic transport arises in a large variety of contexts; these include not 

tG • only astrophysics [T] and cosmology [2J, but also plasma physics [3] and heavy ion 

collisions [3], and even condensed matter physics [51(5] where transport, allbeit non 
relativistic, occurs at bounded speed [5]. There are three main types of models for 
relativistic transport: the purely macroscopic, so-called hydrodynamical models [7], 
the models based kinetic theory [8] , and the stochastic models [9] [10] . The purely 
macroscopic models have been developped since the 1940's [TT], but have serious 
limitations. In particular, traditional Landau-Eckhart theories have been proven to 
\Q , be non-causal |12| and strong arguments |13) . both mathematical and physical, exist 

against the more recent macroscopic theories based on extended thermodynamics 

O; [TU. 

As explained below, the main tool of relativistic kinetic theory is a relativistic 
generalization of Boltzmann famous transport equation. Various relativistic Boltz- 
mann equations have been proposed since 1940. A systematic treatment has been 
proposed only recently in [TSl US] (see also [T7] for a relativistic generalization of 
the Vlasov equation). The treatment proposed in |17| is covariant, but not mani- 
festly so. On the other hand, |15 1 IT6 ] offers several equivalent, manifestly covariant 
^ ' approaches, but fails to offer a purely intrinsic presentation. The relativistic Boltz- 

mann and Boltzmann- Vlasov models have long be assumed |13j to be causa]j, but 
there is, to the best of our knowledge, no formal proof of this assertion in the 
existing litterature. 

Physically realistic stochastic models have been developed since 1997 [18] and 
relativistic stochastic processes now constitute a rapidly expanding field in both 
mathematics and physics. Recent references are e.g. [19l [20] [211 E2] and (23] [24] 
I25[ I26| . where the reader will find intrinsic, covariant and manifestly covariant 
approaches to relativistic diffusions. All relativistic stochastic processes studied so 
far are causal. 
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2 I. BAILLEUL 1 AND F. DEBBASCH^ 

The aim of the present article is threefold. 

(i) propose a clear, intrinsic presentation of the relativistic Boltzmann equa- 
tion, 
(ii) use this intrinsic presentation to establish, for the first time, a clear link be- 
tween the two main branches of relativistic transport models i.e. relativistic 
kinetic theory and relativistic stochastic processes, 
(hi) use that link to offer a simple proof that the relativistic Boltzmann is causal. 

All results are presented in an arbitrary oriented and time-oriented space-time. 

The material is organized as follows. Section 2. offers a presentation of the 
physical aspects of relativistic kinetic theory. Section 3 sets up the geometrical 
tools while Section 3 reviews the definition and properties of the relativistic one- 
particle distribution. Section 4 presents an intrinsic probabilistic interpretation of 
the relativistic Boltzmann equation; the proof that this eqaution is causal is also 
outlined in this section, all technical details being relegated to the Appendix. 

2. Physical aspects of relativistic kinetic theory 



The traditional Boltzmann equation [271 [28] aims at describing the out-of-equilibrium 
dynamics of a dilute gas of non relativistic and non quantum point particles. This 
equation is today best derived from the so-called BBGKY hierarchy [2"51l28| . Sup- 
pose there are N particles in the gaz; according to classical mechanics, the evolution 
of the gaz is then completely determined by a system of 6-/V ordinary differential 
equations fixing the dynamics of the particle positions and velocities. The prob- 
ability of finding, at a certain time t, any k < N particles at a certain point of 
the fc-dimensional phase-space then admits a density with respect to the Lebesgue 
measure in this phase-space and this density is obviously the product of Dirac 
distributions. 

This description of the gaz is of course of little use because N is very large. It is 
also of little physical interest, because the positions and velocities of the individual 
gaz particles cannor be measured. What one observes are rather smoothed out or 
averaged quantities. One ususally supposes that the averaged probability of finding 
at a certain time any k < N particles at a certain point of the fc-dimensional phase 
space, still admits, for all k, a density with respect to the Lebesgue measure in this 
phase-space. It is then straightforward to deduce from the equations of classical 
mechanics, a system of equations obeyed exactly by all these densities. This system 
constitutes what one calls the BBGKY hierarchy. 

If only pair-interactions are taken into account, the equation k of the hierarchy 
is an integro- differential transport equation fixing the dynamics of the fc-particle 
density in terms of the k + 1-particle density. The Boltzmann equation is deduced 
from the k = 1-equation of the BBGKY hierarchy by taking into account the 
high dilution of the gaz and postulating that the two particle density is completely 
determined by the one particle density, and by supposing all interactions between 
particles to be close range interactions which can thus be assimilated to point 
collisions. Note that this restriction does not preclude interactions of the particles 
with an 'exterior' field, independent of the gas dynamics. 

How much of the above picture can one generalize into a relativistic model of 
transport? Unfortunately, not much. Indeed, a consistent relativistic description of 
N interacting particles involves, not only the particle degrees of freedom (positions, 
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velocities), but also the interaction field degrees of freedom. In other words, the 
'mechanical' equations are, in the relativistic regime, a set of differential equations, 
not for 6iV , but for an infinite number of degrees of freedom. Extending the above 
approach to the relativistic realm would thus necessitate introducing densities in 
an infinite dimensional space. Writing a relativistic equivalent of the BBGKY 
hierarchy thus seems a rather formidable task, and dealing with such a generalized 
hierarchy appears even more daunting. 

The classical Boltzmann equation however, taken by itself, can be generalized 
to the relativistic regime; this is possible because the Botzman equation considers 
only one particle densities, which admit a nice relativistic formulation, and treats 
all particle interactions as point collisions. 

The relativistic Boltzman equation is then built in two steps. Step 1: it is possible 
to introduce a natural geometrical object which generalizes to the relativistic realm 
the standard notion of one-particle density. This will be called the relativistic one- 
particle density. Step 2: Consider a gaz of N non quantum but relativistic particles 
immersed in an 'exterior' gravitational field, independent of the dynamics of the 
gaz. Suppose also that all gaz particles interact only through collisions and that 
the gaz particles are free between their collisions i.e. follow geodesies of the exterior 
gravitational field. It is then possible to write down fomally a relativistic Boltzman 
equation obeyed by the relativistic distribution function. The left-hand side of this 
transport equation is simply the action of the geodesic flow on the relativistic one- 
particle density and the left-hand side is a collision term, the general form of which 
is independent of the detail of the particle interactions. 

This article makes clear the probabilistic content of the relativistic Boltzmann 
equation by introducing a well chosen random dynamics and by showing that the 
distribution function obeys Boltzmann equation if, and only if, it is the density of 
the invariant measure of the random process. The probabilistic point of view also 
makes the causal character of the relativistic Boltzmann equation very clear. 

The use of probabilistic methods to investigate the classical homogeneous non- 
relativistic Boltzmann equation is not new and dates back to Kac's suggestion |30) . 
followed by McKean's work [21] . The first breakthrough came from Tanaka's work 
|32) who proved some exponential rate of convergence to equilibrium for Maxwellian 
gases by using probabilistic tools. The enormous industry which followed (|33|. [34], 
[35], [36] etc.) culminated recently with the results by Fournier et al. [37], [38] who 
obtained some uniqueness results for some singular collision kernels by probabilistic 
methods based on couplings. It is however, to the best of our knowledge, the 
first time that a probabilistic view is given on the relativistic Boltzmann equation. 
More generally, it has become more and more clear that the study of random 
processes with values in Lorentzian manifolds can bring interesting insights on 
different questions ranging from the irreversibility problem in relativistic statistical 
mechanics [39] to plasma physics [3D] , [22] , to the study of spacetime singularities 



3. Geometrical setting 

Let (M, g) be a Lorentzian manifold, oriented and time-oriented. Denote by TM 
its tangent bundle, with generic point <p = (m,rh), and by T-'-M the unit future- 
oriented bundle over M. 
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Denote by Volm the volume form on M associated with the Lorentzian metric 
(for which VoLMi(eo, . . . , 63) = 1, if (e<j, . . . , 63) is an orthonormal basis of M at 
some point). Identify the volume form and the volume measure VohM(dm). The 
tangent bundle TM inherits from the Lorentzian structure of M a volume measure 
VohTu(dip) which is the semi-direct product of Volm by the Lebesgue measure 
LEB m (dm) in each fiber T m M, normalized to assign measure 1 to any hypercube 
of T m M constructed on an orthonormal basis: 

VoLTM{d<p) = LEB m (drh) (g> VoLfti(dm), <p = (m, rh). 

At any point m <E M, the metric g m on T m M induces on each hyperboloid T T ^M 
a Riemannian metric; denote by VoL^(dm) its associated volume measure, where 
dm is understood here as a surface element in T^M. The volume measure Vol^im 
is 

VohTi-wiidtp) = Vol^ (dm) (g> VoLu(dm), <p = (m, rh). 

As is well-known, geodesic motion induces a dynamics in TM which leaves the 
bundle T l M. stable: freely falling particles have a velocity of constant norm. Let 
denote by Hq the vector field on TM generating the geodesic motion. Given a 
local coordinate system x : U C M n- R 4 , any tangent vector rh of T m M can 
uniquely be written Y^i—o^^ii w hh the usual notations. The map (m,rh) £ 
TM \-t ((V)ck^3, (w l )o^i^3) defines a local coordinate system on TM. In these 
coordinates, the geodesic vector field Hq reads 

3 



E 



d d 
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dx l %3 dm k 

4=0 

where the L's are the Christoffel symbols of g. 

Given a spacelike hypersurface V, write T X V for {(m,m) G T 1 M; m 6 V}. 
This bundle inherits from the Lorentzian metric a natural volume measure Vol t i v 
which is the semi-direct product of the Riemannian volume measure on V and the 
Riemannian measure in each hyperboloid T^M, for m <G V. Note that TT^V is not 
the unit tangent bundle to (V, g). 

4. One particle distribution function 

Follow the random motion of a typical particle of a relativistic gas, in a spacetime 
(M, 0); it describes a random path ip s — (m s ,m s ) in T X M, wherqj s is a multiple 
of the proper time of the timelike path (m s ). Without loss of generality, we can 
restrict ourselves to the case where m s has unit norm and s is the proper time of 
the particle. The random path can be thought of as a succession of (potentially 
infinitesimal) geodesic segments separated by points where random collisions change 
the velocity of the particle. 

Statistical physics [411 [T51 14"21 [T5] suggests the following assumptions about this 
process. 

• One can associate to any spacelike submanifold V a measure /Li T iy(dy>) on 
the unit tangent bundle T 1 V over V, to be understood as the distribution 
of a typical gas particle hitting V. The measure fijiy(dip) has a density 
fv(<p) with respect to the natural volume measure VoLfiy(d(p) on TV. 



Shocks keep the norm of the velocity constant. 
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• Given any point ip = (m, rh) £ T 1 M, define V v as the set of spacelike 
hypersurfaces V of M containing m and orthogonal to rh at m. The value at 
if of the density /y does not depend on the arbitrary choice of hypersurface 
V £ V v , so fy(<p) is a well-defined scalar; denote it by f(ip) = f(m, rh). 

• At any point m in spacetime, the vector field j(m) = J ■ pT1 M rhf(m, rh) VoL* n (dm) 
represents the particle current at point m. In particular, given any spacelike 
submanifold V with future unit normal toyim) at m £ V, the component 

of j(rn) normal to V at point m represents, for someone who considers V as 
3D space, the 3D or spatial particle density n(m), defined with respect to 
the natural volume measure VoLy (dm) induced by the Lorentzian metric 
of V. Such an observer will also consider T X V as 6D phase-space and will 
thus view g(rh, zay(m))f(m, rh) as the phase-space particle density, defined 
with respect to the natural volume measure on T 1 V. 

The function / is usually called the particle density in T X V, or one particle 
distribution function of the gas. This density completely determines, through 
its first moment j, the particle content of the spacetime. The stress-energy content 
is determined through the second moments of /. Note that the zeroth moment of 
/ has no physical interpretation. 

5. Evolution equation for the one particle distribution function 

Let us associate to the one particule distribution function of the gas a Markov 
process performing geodesic motion in between shock times where it is hit by par- 
ticles of the gas, resulting in a change of its speed. Parametrize this process by 
the proper time of its trajectories in M. The rate at which the shocks happen is 
supposed to depend only on the one particle distribution function and on the cho- 
sen model for the collision mechanism of pairs of particles. Given two particles at 
location m £ M, with velocity rh and rh', denote by p and p' {£ T^M) the outcome 
of the collision of the two particles corresponding to the scattering angle 9 £ § 2 . 
We denote classically by W(m ; rh, rh'; 9) the collision kernel, which represents the 
rate at which the above collision holds; it has the symmetry property 

(1) W(m ; rh, rh'; 6) = W(m ; p,p'; 9), 

for all m £ M and rh, m! £ T m M, usually refered to as the microscopic reversibility. 
See e.g. the book [43j of Cercignani and Kremer for precise models of collision 
mechanisms and collision kernels; these practically important details are irrelevant 
for us in this work. We define our Markov process by its generator: 

(2) 

(Qh)(m,rh) — \Hohj(m,rh) 

{h(m,p) — h(m, m)} W(m ; rh, rh' ; 9) d9 f(m, rh')Voi} m (drh'). 

T ra MxS 2 

Recall that the outcome p of the collision is a function of incoming momentam, rh' 
and the scattering angle 9. Note that as the total rate of collision 

W(m ; rh,rh! ; 9) d9 f(m, rh')Voh m (drh') 

T m MxS 2 
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might be infinite, one would really need the sophisticated tools of stochastic calcu- 
lus to justify the above intuitive picture of the motion as geodesic trajectories in 
between shock times, as these shocks times would not be discrete in case the above 
integral is infinite. This is not our purpose here, though. 

5.1. The relativistic Boltzmann equation. One can associate to any ip = 
(m, m) G T X M, the set V v of spacelike hypersurfaces V of M containing m and 
orthogonal to m at m. Suppose now the hitting distribution by the process of any 
spacelike hypersurface V has a density gv with respect to Vol t i v . The value of 
the density gy at point p will not depend on the arbitrary choice of hypersurface 
V G Vy. Indeed, gy{ip) is the limit of the ratio of the mean number of particles 
hitting a neighbourhood U of ip in TT^V by the volume of that neighbourhood, as it 
decreases to {</?}. Given another V in V v , we can map Wtoa neighbourhood U' of 
<p in T 1 V by a diffeomorphism arbitrarily close to the identity since V and V have 
the same tangent space at m, provided U is small enough. The two limit ratios 
gy(ip) and gy>(<p) will thus have the same value. So gy(<p) = 9(f) i s a natural func- 
tion (scalar) on the unit tangent bundle T'-'-M, named one particle distribution 
function of the Markov process. The function g enjoys the following crucial 
analytical property, proved in Appendix. 

Proposition 1. We have: Q*g = 0. 

So the measure g(p)Vohxim(df) is invariant for the random dynamics, and the 
equation Q*g = is a detailled balance equation. In a more concrete way, proposi- 
tion 1 means that the integral 
(3) 
/ 9(f) \ (Hoh)(ip)+ / {h(m, p)—h(m,rh)}W (m ; rn,m ; 9) d9 /(m, rh')VoL m (drh') > Voh T i M (dip) 

JT 1 M *■ JT m MxS 2 ' 

is null for any smooth functions h with compact support. We write here ip = (m, rh) 
for a generic element p G T 1 ^!. The symmetry property (JlJ of the collision kernel 
and an integration by part^o enable to re-write ((3|) under the form 

(-Hog + C(f, g)) (<p)h(<p) VoL Tm (dp) = 0, 

where 

C{f,g)(f)= / / {g(m,p)f(m,p')-g(m,rh)f(m,m')}W(m-, rh,rh'; d)d6Voh m (drh'), 

that is Hgg = C(f,g). Boltzmann's fundamental chaos hypothesis is equivalent to 
saying that the one particle distribution function of the gas and the one particle 
distribution function of the Markov process coincide: g = j ' . Equation 

H f = C(f,f) 

is the usual form of the relativistic Boltzmann equation. Consult (44] for a totally 
different and axiomatic presentation of the general relativistic Boltzmann equation. 



The vector field Ho has an L 2 (VOL T i M )-dual equal to — Ho as it preserves Liouville measure on 
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5.2. Causal character of the relativistic Boltzmann equation. We show 
in this section how the introduction of the above random dynamics leads to a 
clear understanding of the causal character of the general relativistic Boltzmann 
equation, through proposition 1. We refer the reader to the works [351 1351 R7| and 
|48| of Dudynski and Ekiel-Jezewska for mathematical works on that question in 
the special relativistic case. 

Fix an open spacelike hypersurface V and denote by D + (Y) its future domain 
of dependence: it is the set of points m of M such that any past-directed tiemlike 
path started from m hits V. This set is known to be globally hyperbolic, [49J. The 
next proposition holds for all globally hyperbolic spacetimes although we state it 
for D+(V). 

Proposition 2. One can associate to any point m of D + (V) a positive constant 
T(m) such that any past directed timelike path started from m, parametrized by 
its proper time, hits V before time T(m). 

Proof - It suffices to take for T(m) the length of a future-directed maximal geodesic 
from V to to, whose existence is guaranteed by the global hyperbolicity of 
D + (V) - see e.g. prop. 2.33 in Senovilla's review [50], or consult [5JJ. > 

Consider the T^M-valued Markov process {ip s )s^o = {m Sl m s ) s ^,Q with generator 

g*h = -H h + c(f,h)- > 

it has past-directed timelike paths. Start it from a point (to, to) G T X M with m G 
D + (V). Since / is C/*-harmonic (by proposition 1), the random process (f(ip s )) >n 
is a non-negative martingale. Denote by H the hitting time of T 1 V by (ips)s^o'-, it 
is ahnost-surely bounded above by T(to), by proposition 2. One can thus apply the 
optional stopping theorem and get 

f(m, m) = E( m ,m) [f{i>H )] ■ 

This identity proves the first part of the following statement. Write T 1 D + (V) for 
{(to, to) efM; m G £>+(¥)}. 

Theorem 3. Let (M, g) be any Lorentzian manifold and V be a spacelike hyper- 
surface. The one particle distribution function of a gas is a causal function: its 
values on T 1 D + (V) are determined by its values on T X V. The restriction of / to 
T X V is the minimal set of data needed to determine / on T 1 D + (V). 

Proof - The second part of the statement directly comes from the fact that the 
distribution of tpn has support in the whole of T 1 (l~ (itlq) n V) for a process 
started from the point ipo = (mo, Too)- » 

Appendix 

The result of proposition 1 comes from Kolmogorov's forward equation for the 
transition semi-group of a general Markov process X; we recall it here. 
Denote by x a generic element of the state space of the process and write 

P t (x,h)=E x [h(X t )] 

for the expectation of h(X t ) for a process started from x; write Pt(x,dy) for the 
associated kernel on the state space. Write, as above, Q for the generator of the 
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process. Kolmogorov's forward equation comes from the semi-group property of 
the kernels Pt(x, •), encoded in the Chapman-Kolmogorov equation 

Pt+s (x, h) = / Pt(y,h)P s (x,dy), Vs,£ ^ 0, x in the state space, 

and reads (see e.g. Chap. 1 of [52J) 

(4) j t P t {x,h) = P t {x,Qh). 

In a context where the kernels Pt(x, •) are given by a density p t {x,y) with respect 
to some reference measure dy, equation Q re-writes 

(5) —pt(x,y) = g*p t (x,y), 

where Q* acts on y and is the dual of Q in h 2 (dy). Note however that there is no 
need of densities to make sense of equation (|4]). 

The result of proposition 1 is local in M; it will come as an application of equation 
([5]) by reparametrizing locally the trajectories of the process by a time function 
defined locally on M. The following local construction will be used to that end. 

a) Normal variation of a spacelike hypersurface. Let V be a relatively com- 
pact spacelike hypersurface of M. For m £ V and e € M small enough, define $ e (m) 
as the position at time e of the geodesic started from m, leaving V orthogonally 
in the future direction with a unit speed. Then there exists, as a consequence of 
the local inversion theorem, a positive constant n and an open set U C M such 
that the map $ : ( — 77 , 77 ) xV->M, (£,rn) ^ ®s( m ), is a diffeomorphism. Let us 
further suppose r\ and V small enough for hi to be strongly causal. Writing V e for 
$ e (V), the map $ is the identity on V, and <9 £ $ £ (m) g T^ , m \M is orthogonal to 
T$ e ( m )Y e . The family of spacelike hypersurfaces {^e}ee(-n,n) ^ s ca H e d the normal 
variation of V. The following related notation will be useful. 

Notations. We define a vector field w on U as follows. Given a point m <E Y £ , 
denote by w(m) the future unit timelike vector orthogonal to T m V e ; set w((p) :~ 
w(m), for ip = (m,m). 

• 7 — 7"(v) : ~ ff( tn ( l f 9 )j™) w iM be a function of ip = (m,rh) in the tangent 
bundle of 14. 

• The * TVe -operation will stand for taking the L 2 (VoLTv e )- dual and the *- 
operation for taking the L 2 (VoLri M )-(fuaL 

• For clarity, all objects defined on V or TV will have a hat on them: (p,G,ip £ ... 
defined below. 

• Last, H e will denote the hitting proper time o/TV e C TM. 

b) Reparametrization of the trajectories of the process. Given a point ip 
of T^M take a relatively compact spacelike hypersurface V containing in and do 
the preceding construction. To prove that Q* g = at <p we only need to consider 
what happens near tp. Let us then work in the tangent bundle of Li, where we can 
use the parameter e as a time parameter rather than using the proper time of the 
random trajectories. That is, consider the re-parametrized process {'4>H c }se(--n.ri)] 
it has generator 7 _1 Q. Decompose this operator as follows 

(6) V^ = <M£)eV £ , ^L(<p) = (wf)(<p) + G(f°<Pe)m = (™f)(<p) + {Gf)(<p), 

7 
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where Q is an operator on TV, and where, as a consequence, Q acts only on TV e . 
Now, define the TV- valued process {i> 6 } p , -. := {^e 1 (' l pH e )} (Z( _ s and denote 

by £ E its time-dependent generator. 

c) Proof of proposition 1. Without loss of generality, one can assume that po 
has a smooth density with respect to Voltv and denote by p £ the density of the 
distribution of ij) s with respect to Voltv- By Kolmogorov's forward equation, it 
satisfies the equation 

d £ p e = e^p 6 , 

for all e 6 (—77,77). The operator ^* TV stands here for the L 2 ( VoLTv)-dual of £ e . 
Let us now denote by Vol™ the pull-back on TV by <j> e of the measure Voltv 6 on 
TV £ , and denote by G e its density with respect to Voltv- Then (p e has a density 

G F 



p e = -p^- with respect to VoL^yj it satisfies the equation 



(?) d E fJ. E + -pr- p 6 = V Pe- 

We have written here £* TV ; e g for '- g e9J . Denote by /x e the density of ipH e with 
respect to Voltv e , and consider p and G as functions of e and <p S TV 6 , that 
is, consider them as functions defined on the tangent bundle of U. By its very 
definition, the function p and the one particle distribution function of the process 
are linked through the relation 

(8) n{<p) = lg(v), 

discussed in the third point of section 0J Equation ([7]) can be written in terms of p 

as 

,„s roG — ..TV, 

(9) wp + — -p = Q p. 

G 

The operator Q has been introduced in equation ((6]). It is useful at that stage to 
remark that we have 

as a consequence of the change of variable formula, and since we have a normal 
variation of V. The following lemma is needed to make the final step. 

Lemma. We have for any smooth function / 

Proof - As above, this is consequence of the change of variable formula and the fact 
that we have a normal variation of V. Write T l U for the future unit tangent 
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bundle over hi and take h a smooth function over T X U with compact support. 
(w*f) ((f) h(<p)Voh T i u (d<p) = / f(ip)(mh)(<j))VoL T i M (d<p) 

f(e, (p) (d E h)(e, <p) G e {(p) VoL TV (d^) de 

(-17,17) JT¥ 

(10) =-/ (def){e,<p)h{e,ip)G e ((p) VoL TV (d^)de 

J(-t|,if)xW 

(fh)(e, ip) d e G e {(p) VoL TY (dfi) de 

(-17,17) XTV 

zuf + —r ) (</?) ft(y) VoL T i M (dv3). 

t> 
As a consequence of this lemma we can use the decomposition given in equation 
© to write equation ^ as G*(^) = 0, that is Q* g = 0, by equation (8). 
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